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ABSTRACT 

We compute the fluctuations in gravitational lens image positions and 
time delay caused by large scale structure correlations. We show that these 
fluctuations can be expressed as a simple integral over the density power 
spectrum. Using the COBE normalization we find that positions of objects at 
cosmological distances are expected to deviate from their true positions by a 
few arcminutes. These deflections are not directly observable. The positions 
of the images relative to one another fluctuate by a few percent of the relative 
separation, implying that one does not expect multiple images to be produced 
by large scale structures. Nevertheless, the fluctuations are larger than the 
observational errors on the positions and affect reconstructions of the lens 
potential. The time delay fluctuations have a geometrical and a gravitational 
contribution. Both are much larger than the expected time delay from the 
primary lens, but partially cancel each other. We flnd that large scale structure 
weakly affects the time delay and time delay measurements can be used as a 
probe of the distance scale in the universe. 

Subject headings: gravitational lenses — cosmology: large-scale structure of the 
universe 



1. Introduction 
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The possibility of studying the physical parameters of the distant universe using 
gravitational lenses (GL) was first suggested in the 1960s. In particular, Refsdal (1964, 
1966) suggested that one could determine the masses of galaxies and the Hubble constant 
using the observed image properties, most notably their positions, magnifications and 
time delays between images of the same source. The latter became especially interesting 
after the time delay in the system 0957+561 was measured (e.g. Vanderriest et al. 1989; 
Lehar et al. 1992) and a value of Hq derived (Rhee 1991; Roberts et al. 1991). Alcock & 
Anderson (1985, 1986), Watanabe, Sasaki & Tomita (1992) and Sasaki (1993) criticized the 
method, arguing that large scale structure (LSS) might significantly affect the time delay. 
Unfortunately, their arguments were only qualitative and could not give realistic predictions 
of the amplitude of fiuctuations. For example, Falco, Gorenstein & Shapiro (1991) found 
that the effect Alcock & Anderson discussed can change the derived value of Hq by only 
a few percent. In light of this many workers in the field have taken an optimistic view 
and assumed that the derived value of Hq gives at least an upper limit to the actual value 
(e.g. Borgccst & Refsdal 1984). These arguments are based on the fact that mass density 
is always positive and therefore always focuses the rays. However, this is only correct in 
Newtonian gravity and becomes invalid in cosmological applications, where underdensities 
such as voids give an effective negative mass density (e.g. Nityananda & Ostriker 1984). In 
general, the question whether LSS could significantly affect the measured properties of the 
lens has remained largely unanswered. 

In this paper we present a calculation of position and time delay dispersions in a GL 
system. The calculation is done using linearized general relativity for some reahstic COBE 
normalized cosmological models. In §2 we present the method, based on the geodesic 
equation for a perturbed Robertson- Walker metric. A similar approach has been used by 
hinder (1990) and Cayon, Martinez-Gonzalez & Sanz (1993a, 1993b) to study the GL 
effects on the cosmic microwave background and by Kaiser (1992) to derive the ellipticity 
correlation function of distant galaxies. An alternative method, based on optical scalars, has 
been developed by Gunn (1967) and applied to the ellipticity and magnification correlation 
function calculations by Babul & Lee (1991), Blandford et al. (1991) and Miralda-Escude 
(1991). In §3 we apply the method to compute the fiuctuations in the image position 
relative to the true position and relative to another image position. In §4 we present the 
calculation of fluctuations in time delay. In §5 we present the conclusions and comment on 
their agreement with previous work on this subject. 



2. Formalism 
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Our framework is a perturbed flat Robertson- Walker model with small-amplitude 
scalar metric fluctuations. In the longitudinal gauge (Bardeen 1980; Mukhanov, Feldman 
and Brandenberger 1992) one can write the line element as 



ds = a 



2/ 



-(1 + 20)rfr2 + (1 - 2(p)dx ■ dx 



'V. 



We assume that |0| <^ 1 and neglect all the terms of order 0{(j)'^) and higher. This is a good 
approximation almost everywhere in the universe, except near black holes. We adopt units 
such that c = 1. In the line element above we neglected the contributions from vector and 
tensor modes and the gravitational effects from anisotropic stresses. These approximations 
are valid nearly always, especially in the regime of interest for us, which is the matter 
dominated era with fluctuation wavelengths small compared to the Hubble distance. In this 
case can be interpreted as the Newtonian potential and, neglecting the contributions from 
wavelengths larger than the Hubble distance, it obeys the cosmological Poisson equation 

VV = 47rGa2(p-p), (2) 

where p is the local density and p the mean density in the universe (e.g. Bertschinger 1993). 
We denote the time dependence of the potential with F{t), which is independent of the 
scale in linear perturbation theory, assuming the dominant matter component has negligible 
Jeans length. For the zero curvature model with a cosmological constant, Q^n + = 1) one 
has (Heath 1977) 

^^^^ J^O^a 

where X = a/{^lm + ^\a'^) and HqT = da/{nma + ^lxa'^Y^'^; we normalize a(r) = 1/(1 + 2;) 
to a = 1 today. Hq is the Hubble constant and Qm and Qx are the matter and vacuum 
densities, respectively, in units of the critical density. The simplest model has Qm = 1, for 
which F(r) = 1 and a(r) = {HoT/2y. 

Suppose a photon is emitted from a distant source toward the observer (see Figure 
|I]). If there are no perturbations present the photon will travel along a null geodesic in 
the radial direction ^ , x{r) = rf and r(r) = tq — t, where tq is the conformal time today. 
Adding a perturbation changes the photon trajectory into x{r) = rr + Xj_{r), where 
x±{r) is the excursion of the photon in the direction orthogonal to f. The 3-tangent to 
the geodesic curve is given by n = dx/dl, where dl = (dx ■ dxY^"^ is the comoving path 
length along the geodesic. Null geodesies obey the relation dr = {1 — 2(j))dl. Even when 



^Here and throughout the paper aU 3-vectors are defined in the unperturbed comoving coordinates of 
3-space, which is a hypersurface of constant r. The geometry of this 3-space is a simple Euclidean geometry. 
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metric perturbations are present, we can continue to parametrize the geodesic with the 
unperturbed radial coordinate r. The relation between r and / is given hj dr/dl = n ■ r. 
The linearized space part of the photon geodesic equation derived from the metric in 
equation (|I|) gives the rate of change in photon direction, 

d/Tt 

— = 2n X (n X V0) (4) 
dl 

(e.g. Weinberg 1972 eq. 9.2.6-7; note that there is a factor of 2 missing in eq. 9.2.7). 
For weak gravitational fields (|0| ^ 1), can be viewed as providing a force deflecting 
the photons and affecting their travel time while they propagate through unperturbed 
space-time. We shall adopt this quasi-Newtonian interpretation throughout this paper. 

Because n is normalized, it is sufficient to consider only its two components orthogonal 
to r. We define the two-dimensional photon direction angle 7(r) = [7i(r), 72(r)] relative to 
the unperturbed photon direction — f (taken as the third direction) as n = (71,72, —1). We 
assume small deflection angles (this will be justified in §3). 

The evolution of 7(r) is given by 

j{r) = j{ros) + a_L(r, ros), (5) 

where 'yiros) is the initial photon direction (here vqs is the comoving distance between the 
observer and the source) and ct±{r,ros) = [(^±,i{r,ros),C(±,2{r,ros)] is a two-dimensional 
deflection angle produced by the potential along the geodesic between the source and the 
point under consideration at r. From equation (^ follows 

r^os 

OL±{r,ros) = j g±ir)dr. (6) 

We introduced g±{r) defined as 

9±{r) = -2V±<j)[x{r), r(r)] = -2[V - f (V • fMx{r), r(r)]. (7) 
The photon excursion in the plane perpendicular to f is given by 

pros 

x±{r) = / j{r)dr. (8) 

Jr 

The initial photon direction "flros) must be chosen so that a;j_(0) = 0, i.e., the ray must 
pass through the observer's position. This gives the lens equation 



/ 'y{r)dr = 
Jo 
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or 



pros 



roslijos) = - OL^{r,ros)dr. 



(9) 



This lens equation is valid for an arbitrary mass distribution between the source and the 
observer. It cannot be solved explicitly in general, because Q:j_(r, ros) depends on the 
unknown initial photon direction 'y{ros)- Instead, one has to solve an integral equation 
using, for example, the ray-shooting method. 

We may also ask what is the photon direction at the observer's position so that 
x±{r-os) = 0. This is given by the image position angles 7(0). For the special case of a 
single thin lens one has 



where 6l is the bending angle in the lens plane and tol is the comoving distance between 
the observer and the lens. Equation (|^) then gives 



which is the usual lens equation in the thin lens approximation (e.g. Schneider, Ehlers & 
Falco 1992; Blandford &: Kochanek 1987). Here r^s is the comoving distance between the 
observer and the source and tqs = fOL + r^s- 

The distances can be calculated for a given cosmological model from the source redshift 
Zs and lens redshift zi (assuming that the two redshifts can be measured). For example, for 



= lone has roL = 2H^^[l- {1 + zi)-^'^] and r^^ = 2H^\l + zi)-^'^ - {1 + Zs)-^/\ We 



use the unperturbed comoving distance-redshift relation, because the deviations from it are 
0(0) and can be neglected to the first order. Note that we expressed the lens equation using 
comoving distances. One can reexpress it, if one wishes, with angular diameter distances 
using the relation (Ias = ^As/iX + ^5)5 where zs is the redshift of the source and (Ias and 
tas are, respectively, the angular diameter distance and comoving distance between the 
point A and the source. 

Similarly, the net time delay along the photon path relative to the unperturbed path is 
given by 




(10) 



^os7(0) 



(11) 




(12) 



where 7^(r) = 7i(r) + 7|(r) and we used the small angle approximation; the unperturbed 
null geodesic equation may be used for r(r) = rg — r. The first and second term in equation 
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(p^) are the geometrical and gravitational (also called potential or Shapiro) time delay 
contribution, respectively. Again, for the particular case of the thin lens approximation this 
reproduces the usual time delay equation, as shown in §4. 

Although we have derived our equations using the geodesic equations of a perturbed 
Robert son- Walker metric, the final expressions agree with the usual thin lens equations, 
when we interpret as the Newtonian potential, obeying equation ^ In addition, one does 
not have to assume that density perturbations are small, as long as \(f)\ -C 1. These are 
advantages of the longitudinal gauge compared to other gauge choices (e.g. synchronous 
gauge). The particular strength of the approach presented here is that it remains valid 
in cosmological apphcations and can give the deflection angle and time delay due to any 
matter distribution between the source and observer. Use of comoving coordinates and 
conformal time greatly simplifies the GL equations. In particular, all the equations above 
use a simple Euclidean spatial geometry and there is no need to use angular distances or 
worry about expansion factors. The general expressions above would be more cumbersome 
when expressed with angular diameter distances. 

We are interested in LSS effects on the image properties of lenses. Because of the 
statistical nature of cosmological theories one can only predict the ensemble averages of a 
given quantity, such as its mean and variance. All the GL effects are given through the 
gravitational potential 0[a;,r(r)]. Its statistical properties can best be described with the 
Fourier transformed potential 0(fc), 

(t)[x, r(r)] = j d^k(j){k)e'^ ' '^F[r(r)]. (13) 

The ensemble mean and variance of the Fourier transform of the potential are {(f){k)) = 
and {(j){k)(j)*{k')) = P^{k)5'^{k — k'), where P(j){k) is the power spectrum of the potential 
(see e.g. Bertschinger 1992). We will use the power spectrum of the potential in this paper 
because it leads to the simplest expressions, but note that one can easily reexpress the 
results with the density power spectrum using equation 0. 

A particularly convenient normalization of P(p{k) is given by the cosmic microwave 
background anisotropy measured by COBE- the quadrupole Q2 = (AT/T)2 = (6±1) x 10~^ 
(Smoot et al. 1992; Seljak & Bertschinger 1993). On the large scales probed by COBE 
the dominant contribution to AT/T is given by the Sachs-Wolfe (1967) effect, induced by 
the same metric fluctuations that cause the fluctuations in time delay and image positions. 
Assuming no tensor mode contribution to CMB anisotropies and isentropic (adiabatic) 
fluctuations we can express the quadrupole in terms of the power spectrum of the potential 
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(Bond & Efsthathiou 1987) as 

Q2 = -V^/ k'P^{k)jl{2k/Ho)dk. (14) 
9 JO 

Here j2{x) is the spherical Bessel function of order 2 and K2 > 1 is the amphfication 
coefficient due to the time dependent potential (Kofman & Starobinsky 1985). If 
flm = 1 the potential is time independent and = 1. For the scale-invariant 
Peebles-Harrison-Zel'dovich spectrum P(j,{k) = Ak~^ one gets 

Ql = (15) 



3. Fluctuations in Angular Position 



The first question we will ask is what is the rms fluctuation in the photon direction at 
the observer's position relative to the unperturbed direction. This is defined as 

= -^{a±{0,ros) ■ ai(0,ros))'/', (16) 

where ct±{0, ros) is the total defiection angle accumulated between the source and the 
observer, as defined in equation (|^). In order to compute the variance we will make 
an additional assumption that the statistical properties of the potential sampled by the 
perturbed photon geodesic are approximately equal to those along the unperturbed one, 

{(j)[x,T{r)\(j)[x,T{r)\) ^ {(j)[rf,T{r)\(j)[rf,T{r)\). (17) 

This assumption will be discussed later. 

We want to compute statistical properties of the transverse gradient of the potential 
field integrated along line of sight in the direction 7 and weighted with a function g(r), 

pro 

p{l) = -2 / [V^0(r,7)]g(r)dr. (18) 
Jo 

We will use the Fourier space analog of Limber's equation (e.g. Kaiser 1992) to calculate 
the correlation function from the power spectrum, 

rco pro 

C{l=hi-J2\) = {phi)-Ph2))=8n' k'dk P4k,T{r)]q\r)Mkn)dr, (19) 

Jo Jo 
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where Jo{x) is the Bessel function of order 0. Equation 19 assumes that the scales 
contributing to the correlation function are much smaller than the photon pathlength Tq. 
As shown below, this is usually satisfied in most models of LSS for sources at cosmological 
distances. No assumption on the power spectrum has been made and equation can be 
used both in linear and in non-linear regime. To compute the variance we use the above 
expression with 7 = 0, q{r) = 1 and tq = rps- Assuming linear evolution and = 1 we 
obtain 

/ roo \ 1/2 

= iSir^ros P^{k)Pdk] . (20) 

The rms fluctuation in photon direction at the distance ros from the source, a^, is 
related to the rms angular fluctuation of the true source position relative to the observed 
position. To compute this we have to back-propagate all the photons with a fixed final 
direction and ask what are their angular excursions in the source plane. 



1 r^os 



ros JO 



q:_l(0, r)dr 



ros 



r^os 

/ (y.±(0, r)dr 
Jo 



1/2 



(21) 



Integrating by parts the two terms in equation (^) we find that erg is given by a similar 
expression as with g(r) = (1 — r /ros)- Using equation ^ gives 

To get some intuition about the scaling of the amplitude with the parameters we 
will present an estimate of the fiuctuations for a particularly simple power spectrum 
approximating infiationary models with a physical transfer function: 



pm = { .1-7.4 'i:. (23) 




The spectral indices have been chosen to agree with the cold dark matter model in the 
limits of small and large k. Applying this power spectrum to equations ( pO] ) and ( p2D we 
find 

a, = ^ 7Q2{k,rosf'^. (24) 

This result has a simple physical interpretation. For power spectra like in equation (p^) 
the dominant contribution to gravitational defiection of light comes from the scales near 
the turnover position k^^. A photon travelling through a coherent structure of size k^^ will 
be defiected by ^7 ~ 20 ^ where the last relation assumes cf) is scale invariant for 
k < ko and is therefore fixed by the Sachs- Wolfe effect on the Hubble distance scale. Each 
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region of size kQ makes an independent contribution to the deflection. Since tlie individual 
contributions are random, the photon exhibits a random walk with ~ N^^'^Sj, where 
N = k^ros- Numerical factors aside this agrees with equation (p^. A reasonable value for 
the turnover position in the power spectrum is given by k^^ = 10 Mpc. Taking rps = 1 
Gpc as a typical source distance we find cxe ^ 3 x 10"'^(A;o'"os/100)^/^ ^ l'(A;oros/100)^/^. 
This is small compared to 1, which justifies the small deflection angle assumption. It also 
verifies that the pathlengths are not significantly lengthened by the perturbations. 

However, we conclude that the fluctuations in the angular position ag are of the order 
of arcminutes (see also hinder 1990), which means that the true positions of distant objects 
in the universe, such as quasars, differ from the measured positions on average (rms) by 
this amount. These fluctuations arise already because of the linear structures (voids and 
superclusters) and are present even when there are no nonlinear objects (like galaxies and 
clusters) intersecting the photon trajectories. The fluctuations are much larger than a 
typical image separation in a lens system, which is of the order of a few arcseconds. The 
large total deflections are not directly observable in a single lens system, because only the 
relative positions of images can be measured. Figure |^ shows photon propagation in a 
typical two image GL system. Despite the fact that the deflection of any single photon ray 
can be large, the lens equation will still give the same solution as in the unperturbed case, 
provided that LSS deflects the two photons approximately by the same amount. This will 
be examined below. Although the deflection of a photon ray relative to the unperturbed 
direction is not directly observable from the positions, one might worry that it could 
produce significant time delay fluctuations. We will address this question in §4. 

Let us now calculate LSS effects on the relative image positions, by calculating a^-y and 
cr^g, the dispersion in the relative direction and in the relative angular position separation 
between two image rays. We will denote the two rays with A and B, separated in direction 
at the observer's position by an angle A70 = 75^(0) — 7f (0), where the orientation of the 
coordinate system was chosen so that the image separation vector lies in the direction ei 
(so that 7^(0) = 7^(0)). We divide the potential into a stochastic part, which describes the 
LSS, and a non-stochastic part describing the primary lensing object. We assume there are 
no correlations between these two parts. Taking the expectation value of equation (P) the 
stochastic contributions average to and we obtain the usual lens equation in the thin lens 
approximation. 

The difference between the two direction vectors caused by LSS between the lens and 
the observer is given by 

A7(roL) = A7oei + ^^^(O, tol) - ^^''(O, tol), (25) 

where cxj_^{0,roL) and ck_l'^(0, t^l) are the LSS caused deflections for the rays A and B, 



respectively. We excluded the non-stochastic deflection from the primary lens itself. For 
a flxed photon separation angle at the observer's position A70, the rms fluctuation in the 
angle between the two rays at the lens position is given by 



-| 2\ 1/2 



C(A7o) - C7(0) 



nl/2 



(26) 



where C(7) is given by equation (|T9D with q{r) = 1 and tq = roi- Assuming in addition 
that the scales contributing to the fluctuations are much larger than the typical separation 
between the two rays (one can always verify this assumption by evaluating a^-y from 
equation |19]), we obtain 



2vrA7o 



OL 



1/2 



(27) 



Replacing A70 with {roL/fLs)^lo and tol with r^s in equation ^ gives the rms fluctuation 
between the two ray directions accumulated between the source and the lens. Summing 
the two contributions in quadrature gives the rms fluctuation accumulated between the 
source and the observer, neglecting the correlations between paths on either side of the 
lens. Similarly, the dispersion of the angular position in the lens plane is given by setting 
q{r) = (1 — r/roi), which gives for the contribution between the observer and lens. 



10 



27rA7o 



'OL 

30 Jo 



P^{k)k^dk 



1/2 



(2^ 



For the simple power spectrum of equation ( pSf ) we flnd 

= yiOcTAe ~ 4Q2(fcoroL)'/'A7o ^ 0.025(fcoroL/100)='/2^7o. (29) 

Again, there is a simple physical explanation of this result. Two photons separated by an 
angle A70 sample different potentials, Scj) ~ (V^0)r A70 ~ /i;o0^A7o, including only the peak 
power contribution around k^. The separation r7o is largest at the lens, but only falls to 
one-half for distances half and three-halves as far, so it is a reasonable approximation to flx 
r to Vol- A coherent structure of size k^^ leads to an angular difference of (5A7 ^ 250 and 
there are = k^roL random and independent contributions. The total angular difference 
is just an incoherent sum of individual contributions, a^-y ~ A'y, which, numerical 

factors aside, reproduces equation (PUj). 

A more direct way to estimate the amplitude of image position fluctuation is to use 
the observations of correlated distortions of distant galaxy images. This can be described 
by the ellipticity correlation function 6*^^(7) (Blandford et al. 1991), which describes the 
correlations in the ellipticities of galaxy images as a function of angular separation 7. The 
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ellipticity correlation function at zero lag, Cpp{0), can be related to the power spectrum 
using the expression 



where we assumed for simplicity that all the galaxies lie at the same distance Vg (Kaiser 
1992; Blandford et al. 1991). From this one sees that 



For most cosmological models and distances, the linear theory prediction of equation 



effects, which are somewhat uncertain. Theoretical estimates (Kaiser 1992) and N-body 
simulations (Blandford et al. 1991; Miralda-Escude 1991) suggest that Cpp(O) is unlikely 
to exceed 10"'^. This is also supported by the observational data. Mould et al. (1994) 
report a detection of a sig nal with a value Cpp{l' < 7 < 5') = (5.6 ± 0.6) x 10-^ which, 
after seeing correction, implies average ellipticity within a few arcminutes radius of about 
0.05. Assuming this value we find (TA7/A70 ~ 0.08(roL/'rg)^/^, with ?a 0.6Hq^. However, 
the authors could not exclude the possibility that the observed signal is due to systematic 
effects. This is also suggested by Fahlman et al. (1994), who report a null detection of 
average ellipticity within a 2.76' radius aperture with a sensitivity of about 1.3%, which 
after adjustment to Vg above implies an upper limit (TA7/A70 < O-OS^roL/rg)^^"^- These 
measurements give average ellipticities in typically arcminute size windows and do not probe 
Cpp{'~f) on scales below 1'. Observationally it is difficult to give reliable estimates on smaller 
scales because one has to distinguish between the signal and the noise from the intrinsic 
ellipticities of galaxies. A rather weak upper limit on Cpp(O) can be obtained simply from 
the average ellipticity of galaxies, which is of the order of 0.4 and is dominated by intrinsic 
ellipticities. Despite some uncertainty from the model predictions and observations, it 
appears unlikely that the relative fluctuations in the image separation angle exceed a level 
of a few percent. 

The conclusion above justifies the assumption that the rms fluctuation is small 
compared to the measured image separation. The fact that crA7/A7o <C 1 also implies that 
one cannot have multiple images produced by LSS alone. Therefore, multiple images can 
only be formed from nonlinear structures, such as galaxies or clusters. This conclusion has 
previously been obtained using N-body simulations by Jaroszynski et al. (1991) and using 
semi-analytical methods by Bartelmann & Schneider (1991). The fluctuations in angular 
image separation, although small, are in most cases larger than the observational errors 
on the image positions (typically less than 0.01 arcsecond/arcsecond ~ 10~^). Therefore, 




(30) 




(31) 



(p7|) gives (JA7/A70 of the order of a few percent. This has to be corrected for the nonlinear 
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LSS effects are a major source of uncertainty in the true image positions. Observational 
efforts in trying to determine the image position with precisions below 0.01 arcsecond are 
redundant and do not improve the lens reconstruction. This effect fundamentally limits our 
ability to reconstruct the lens potential using the image positions and should be included in 
the modelling of lens parameters, once the actual amplitude of fluctuations is determined 
from the ellipticity correlation function measurements. 

We should also justify the assumptions made in our calculations. One is that the 
statistical properties of the potential along the perturbed path are well approximated by 
those along the unperturbed path, as expressed in equation (|T7|). Taylor expansion of the 
potential gives 

(f){x) ^ (j){rr) + x± ■ V_L0(rf). (32) 

Inserting this into the left-hand side of equation ([T7|), we find that for Gaussian random 
fields the relative correction to the right hand side of equation (|T^ is given approximately 
by {x_Lko)^ ~ [3(52(fco'^OL)^^^]^ ~ 10~^(/co?^ol/100)^. Therefore, the approximation in 
equation (|T^ leads to negligible errors. Another approximation we used was to neglect 
the correlations between LSS and the primary lens. This is justified because the two are 
correlated only over a correlation length distance, which is much smaller than the typical 
pathlength. While there are N uncorrelated regions of size along the photon path, 
only one of those is strongly correlated with the primary lens. The contribution from that 
region can be regarded as being part of the primary lens itself. The error due to this 
approximation is therefore of the order of A^~^ = {koros)~^ 1- 



4. Fluctuations in time delay 



In this section we compute the dispersion in time delay between two images. For this 
purpose it is useful to define the time delay relative to the normal of the lens plane and not 
relative to the source-observer line. We define the lens plane at the lens redshift zi to be 
orthogonal to what would be the source-observer line in the absence of LSS effects (dotted 
lines on Figure Relative to the lens plane normal, the incoming and outgoing photon 
direction vectors in the lens plane are -y^'*"', ^^>o"* and -y^''*"', ^B,out ^^le images A and 
B, respectively. The difference between the incoming and outgoing photon direction gives 
the deflection angles in the lens plane, ol^ and . These can be obtained by modelling the 
lens potential using various observational constraints, such as image magnifications, velocity 
dispersion of the lensing galaxy and/or cluster, positions of other images or arcs, etc. 
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In the absence of LSS the time delay between two images is given from equation by 

= '-^[{itn'-{{iin']-'^{r-^^)- (33) 

Here ip"^, ip^ are the integrals of the primary lens potential for the two rays, 

rroL+<^ 

^= (j){r)dr (34) 

with e/roL ^ 1- Equation ^ is the usual time delay expression in the thin lens 
approximation (e.g. Blandford & Narayan 1986; Blandford & Kochanek 1987; Schneider, 
Ehlers & Falco 1992). Note that the difference between the two outgoing photon directions 
gives the observed image splitting. 

Adding LSS moves both the source and the observer, for a fixed lens plane (Figure 
The time delay between the two rays is now given by 

= \ n {[^^'" - (^OL, r)]2 - [t^'^'^ - a?(roL, r)]^ dr 

+ \ fj"^ {[7^'°"* + «l(r, vol)? - [7^'°"* + (r, tol)?] dr 

- 2 12 [0^(r) -0^(r)]rfr + V^^-7/;^ + y^ [0^(r) - 0^(r)]t;r| . (35) 

This equation is similar to equation (^2|), except that here the geometrical contribution 
is measured relative to the normal of the lens plane. The first two lines in equation ( P3| ) 
give the geometrical time delay between the lens and the source and between the lens and 
the observer, respectively. In the third line we have written the gravitational time delay 
contribution coming from the potential between the lens and the source, from the primary 
lens potential and from the potential between the lens and the observer, respectively. The 
LSS (q:_l, 0) and primary lens (7*", 7°"*, ip) contributions are thus explicitly separated. 

Let us calculate the time delay contribution between the ray A and the fiducial 
ray accumulated between the lens and the observer. The fiducial ray is defined to start 
perpendicular to the lens plane and end at the observer's positions. The total time delay is 
obtained by adding a similar contribution from the lens to the source and subtracting the 
same terms for the ray B. The fiducial ray direction is given by 

T{r) = OL^{r,roL)= 9±{r)dr, (36) 



where g±{r) is computed along the fiducial ray. The direction of the ray A is 

^^(r) =7/(r)+7^'-* + Aa^(r,roL). (37) 
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Here Acx^{r,roL) is the difference between the LSS induced ray deffections at r and can be 
calculated using the Taylor expansion of g± around the fiducial ray, as we did in §3. The 
image ray position relative to the fiducial ray is 

xi{r) = xiiroL) + Y^''^{roL - r) + Aal(r', roL)dr' . (38) 

Jr 

The initial lens plane position of the image ray relative to the fiducial ray, xj_{roL)) is not 
a free parameter, since it has to satisfy the constraint 

xj_{Qi) = 0. (39) 

From this we obtain 



xj_{r) = -7°"*'^r 



r Aai{r',roL)dr'. (40) 
Jo 



The gravitational time delay contribution is obtained from the Taylor expansion of the 
potential around the fiducial ray, 

Atgrav = {g±ir) ■ xiir) + 0[ixiirf]}dr. (41) 

The geometrical time delay contribution is given by 

/ out,A\2 r r 

^ + / ^ ^/(r) . Y-^'^ + y (r) . Aal(r, tol) 
Jo L 

dr. (42) 



2 



2 

Integrating by parts the terms involving 'y^{r) we find 

ri'OL 

/ (r) • [7°"*'^ + Aal(r, roL)]dr = 

Jo 

fj' [oAr) ■ /j7°"*'^ + Aal(r', tolW^ dr. (43) 

This is exactly cancelled by the first order term in Atgrav (equation ^). Therefore, 'j^ 
completely drops out of the time delay expression and using time delay measurements we 
cannot infer any information on the absolute deflection angle. This is quite remarkable, 
given that separately the geometrical and gravitational LSS induced fluctuations are 
approximately 15 yr (fcoro5/100)^/^(ros/lGpc)(A7o/l"), much larger than the expected 
time delay from the primary lens itself, of the order of 0.1 yr (ro5/lGpc)(A7o/l")^. 
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Adding the geometrical and gravitational time delay contributions we finally obtain 

f-,out.A\2 rroL 1 

^^A ^ a ^ / r o«t,A . ^^A(^) ^ -AaiirY + 0{{xi{ry)]dr. (44) 

2 Jo 2 

The first term in equation (H3|) gives the largest contribution and is the term that one would 



also have in the absence of LSS (compare with equation |33|). The second term is smaller 
than the first term approximately by Aa^(0, roL)/7°"* ~ (JA7/A70. The last two terms in 
equation (|^) are further suppressed by crA7/A7o relative to the second term and can be 
neglected. 

What, then, is the LSS induced fluctuation that causes the reconstructed time delay 
to differ from the true time delay? The observer measures the image separation angle 
that is almost, but not exactly, given by 'y^'""* — ^B,out^ \\^q reconstructed time 



delay differs somewhat from equation (|33D . It is not possible to give an exact prediction 
of the reconstructed time delay without specifying the detailed lens model and taking into 
account all of the observational constraints. It is clear, however, that the fluctuation in 
the reconstructed time delay is due only to the fluctuation in the relative angle separation, 
part of which is described by the second term in equation (|^). Given that (TA7/A70 -C 1, 
the relative effects on the time delay will also be of that order. We conclude that the time 
delay fluctuation induced by LSS is of the order of crA7/A7o, which is of the order of a few 
percent for sources and lenses at cosmological distances. 



5. Conclusions 



We investigated the LSS effects on measurable properties of gravitational lens systems, 
in particular on the image positions and time delays. The method we use is based on a 
geodesic equation in a weakly perturbed flat Robertson- Walker metric and is valid for a 
general matter distribution between the source and the observer. The advantage of this 
approach compared to previous work on this subject is that it only assumes the knowledge 
of evolution of density power spectrum, which can easily be related to other measurements 
of LSS to obtain quantitative predictions. The same approach can also be used to calculate 
light propagation in non-flat universes and the conclusions in this paper do not signiflcantly 
depend on the assumed value of Vt. We flnd that the rms fluctuation in the relative positions 
of images is (JA7/A70 ~ 0.025(A;oT'ol/100)^/^, for a LSS density power spectrum peaking at 
wavelength k^. For most realistic models of LSS this is much smaller than unity and so one 
does not expect multiple images generated from LSS. Nevertheless, these fluctuations are 
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likely to be larger than the observational errors and should be included in the modelling 
of lens parameters. Similarly, rms fliictuation in time delay due to LSS are caused only 
by the uncertainties in the relative image positions and are also approximately given by 
0.025(A;o?"o/100)^/^. Therefore, LSS does not significantly affect the time delays. While 
the same method can be used to predict fluctuations in the relative image magnification 
and orientation, a simple estimate shows that the effect on these observables is negligible. 
The rms fluctuation in relative magnification between two images AM/M is given 
approximately by AM/M ~ (V ±M)ros^lo/M ~ (koros^7o)Cpp{Oy/'^ ~ 10~^, well below 
the measurement errors. 

The conclusion regarding LSS effects on time delay measurements disagrees with 
the conclusions reached by Alcock & Anderson (1985, 1986), Watanabe et al. (1992) 
and Sasaki (1993). Sasaki (1993) suggests that the coupling between the absolute and 
relative deflection, which in our language is proportional to a^aA-y, generates large time 
delay fluctuations. As we have shown in §4, this term actually vanishes once both the 
gravitational and geometrical time delay contributions are included. Alcock & Anderson 
(1985, 1986), Watanabe et al. (1992) and Sasaki (1993) have argued that although the 
universe is homogeneous on large scales, individual photons may travel through a region 
where the density differs from the average density in the universe. If one assumes that 
the density in the beam is some constant fraction of the average density, then one can use 
the Dyer-Roeder angular distance redshift relation (see Sasaki 1993 for a comprehensive 
discussion) to deduce the value of the Hubble constant. The result depends on the unknown 
density in the beam and one may conclude that this signiflcantly affects the Hubble constant 
determination from the time delay measurements. For example, Alcock & Anderson (1985, 
1986) argue that because of LSS correlations there may be a significant overdensity in the 
particular direction of the lens. This conclusion neglects the fact that LSS and the primary 
lens are correlated over only a correlation length distance, which is typically much smaller 
than the photon travel distance. Similarly, Watanabe et al. (1992) and Sasaki (1993) 
argue that we might live in a universe where part or all of the mass is concentrated in 
small clumps away from the photon beam, which are unable to significantly affect the light 
propagation. Both descriptions of LSS are oversimplified in that they do not account for the 
stochastic nature of LSS. Dynamical measurements over the last decade show large amounts 
of dark matter present on all scales and so one cannot assume that photons are travelling 
through a uniformly filled or empty beam when away from the primary lens. Moreover, 
the average redshift-angular distance relation should not differ from the homogeneous case, 
once we average over all lines of sight (including those through the clumps of matter). 
If redshift-angular distance relation for most lines of sight is to differ significantly from 
that in the homogeneous universe, then this requires large rms fiuctuation around the 
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mean. This would imply large o^j^^ but, as we argued in §3, no such large fluctuations 
are either expected or observed. Therefore, the use of angular distances as given by the 
homogeneous universe model appears to be adequate. Another reason for the disagreement 
with Watanabe et al. (1992) is that these authors assumed a large contribution to the 
fluctuations from small (galactic) scales. As wc argued, the data at present show little 
support for such a large contribution on those scales, although for a more quantitative 
prediction better measurements of power spectrum on small (subarcminute) scales would 
be needed. 

If LSS does not induce significant time delay fiuctuations, then this would remove 
one of the major objections against using time delay measurements to determine Hq. 
Signiflcant problems related to the robustness of the lens reconstructions still remain 
and are preventing the method at present from giving a reliable estimate of Hq (see e.g. 
Bernstein, Tyson & Kochanek 1993 for a discussion of lens reconstruction in 0957+561). 
Moreover, the above analysis does not exclude the possibility that a homogeneous sheet of 
matter is present in the lens plane, because in our treatment this sheet of matter is part of 
the primary lens. Uniform matter distribution, which is likely to be ovcrdcnse close to the 
primary lens, cannot be determined from the image positions, but it does affect the length 
scale and makes the deduced value of Hubble constant larger than the true value (Borgeest 
& Refsdal 1984; Falco, Gorenstein & Shapiro 1991). Although severe, these problems are 
not unsolvable and GL time delay method remains one of the few methods that can provide 
information on the global distance scale and geometry of the universe. 
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Fig. 1. — Photon propagation relative to the source-observer hne. 



Fig. 2. — Schematic diagram of a typical lensing discussed in the text. Sohd 

lines represent true photon trajectories, dashed lines apparent trajectories as seen from the 
observer's position and dotted lines the unperturbed trajectories as seen from the lens plane 
in the absence of LSS effects. The apparent image and lens positions are denoted by A', B' 
and L', respectively, and can be far from the true positions. 



